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Abstract 

During the last decade, self-affine geometrical properties of many growing 
aggregates, originated in a wide variety of processes, have been well charac- 
terized. However, little progress has been achieved in the search of a unified 
description of the underlying dynamics. Extensive numerical evidence has 
been given showing that the bulk of aggregates formed upon ballistic aggre- 
gation and random deposition with surface relaxation processes can be broken 
down into a set of infinite scale invariant structures called " trees" . These two 
types of aggregates have been selected because it has been established that 
they belong to different universality classes: those of Kardar-Parisi-Zhang 
and Edward- Wilkinson, respectively. Exponents describing the spatial and 
temporal scale invariance of the trees can be related to the classical expo- 
nents describing the self-affine nature of the growing interface. Furthermore, 
those exponents allows us to distinguish either the compact or non-compact 
nature of the growing trees. Therefore, the measurement of the statistic of 
the process of growing trees may become a useful experimental technique for 
the evaluation of the self-affine properties of some aggregates. 

PACS numbers: 64.60.Ht, 05.40.+j, 61.43.Hv, 05.70.Ln 



I. Introduction. 

The study of growth processes has recently attracted considerable attention as a conse- 
quence of many factors such as the formulation of the theory of fractals since the pioneer- 
ing work of Mandelbrot^, the development of new experimental techniques (e.g. scanning 
tunneling microscopy) and the availability of computer facilities for graphical and numeri- 
cal simulations. So far, it is well established that spatial scaling structures originated from 
growth processes are extremely common in nature. In fact, a great variety of systems exhibit 
self- affinities over extended ranges of spatial and temporal scales^^^. Therefore, extensive 
theoretical and experimental research has focused on the characterization of the self-affine 
nature of the structures resulting from growth processes^"^. 

The phenomenological scaling approach to the dynamic evolution of a self-affine interface 
early developed by Family and Vicsek^'^ has become a useful tool to characterize self-affine 
roughness. Considering a fiat, d — 1 dimensional surface at time t — and pointing the 
attention to the growing process that occurs essentially parallel to the surface, it is possible 
to assume without loss of generality that there exists a well defined growth direction and 
that the interface can be described by a function /i(x, t) that gives the height of the interface 
at time t and position x. Of course, such height is measured from the initial flat surface 
at t=0. If the interface cannot be described by a single valued function of x, the function 
/i(x, t) gives the maximum height of the interface at x. Considering a section of the sample 
having a typical size L (in each of the d—1 dimensions of the surface) the average height of 
the interface at time t is defined as: 

where the summation runs over all x's. The interface width w{L, t) at time t may be defined 
by the rms of the height fiuctuations given by: 

w{L,t)^{{h'{t))-{h{t)ff"- (2) 
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The Family- Viscek scaling approach assumes that: 

w(L,t) ^ L'^Fit/L'), (3) 

where F{y) oc for y << 1 and F{y) — > constant for y >> 1, with z — a/ (5. Also a, (5 
and z are the roughness, growing and dynamic exponents, respectively 

The dynamic exponent z describes the evolution of a correlated region with time: initially 
different parts of the interface are independent, but regions of correlated roughness form over 
time and their size grows as ^ oc t^^^. Thus, for a finite sample of side L and t ^ oo the 
width of the growing interface reaches a statistically stationary state so that w{L) oc L". 
Furthermore, the overall width of the interface grows as t^ until it saturates at L". 

In contrast to the progress achieved in the characterization of the self-affine behavior 
of interfaces, little attention has been drawn to the description of the internal structure 
of the growing system. In this work, extensive numerical evidence is presented showing 
that the bulk of two archetypical growth models, namely ballistic deposition and random 
deposition with surface relaxation, can effectively be rationalized on the basis of a treeing 
process; i.e, any growing structure can be thought as the superposition of individual trees. 
These individual trees can be defined as follows: in the case of the ballistic deposition model 
(figure 1(a)), a newly deposited particle is assumed to belong to the same tree as that of the 
nearest neighbor particle where it is attached. If the deposited particle has more than one 
nearest neighbor belonging to different trees such as particle A in figure 1(a), one of them 
is selected at random and the particle is incorporated into that tree . In the case of the 
random deposition with surface relaxation model (figure 1(b)), a newly deposited particle 
is assumed to belong to the tree corresponding to the impingement site, before its eventual 
relaxation. See also figure 1 caption for further details. 

Those trees that spread out incorporating additional growing centers; e.g. capturing 
particles, developing new branches, etc., are said to be "alive" (for example trees 1,2 and 
4 in figures 1(a) and 1(b)). In contrast, some other trees may stop their growth due to 
shadowing by surrounding growing trees, so they become "dead trees" such as tree 3 in 
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figures 1(a) and 1(b). The structure of dead trees is frozen in the sense that it cannot 
be modified by any further growth. In order to determine the distribution of dead trees 
one starts the growing process following the specified set of rules and after some fixed time 
all growing points, i.e. all sites of the aggregate where further growth is still possible, are 
identified as seeds for these trees. Therefore, the number of such seeds becomes the initial 
number of trees. During the subsequent process the competition between growing trees 
dominates the dynamics. Some trees become dead (frozen) and eventually only a single tree 
may remain. 

Within this context, the aim of this work is to perform an extensive numerical investiga- 
tion of the dynamics of evolving trees of two growing processes, namely Random Deposition 
with Surface Relaxation (RDSR) and BaUistic Deposition (BD). Furthermore, tests of the 
scaling relationships relating critical exponents of the treeing process to those describing 
the self-affine nature of the aggregate are performed in dimensions 1 < d < 5. The RDSR 
and BD processes have been selected because they belong to different universality classes, 
namely the Edward- Wilkinson (EW)^'* and Kardar-Parisi-Zhang (KPZ)"^'^, respectively. 

The manuscript is organized as follows: in section II detailed definitions and a discussion 
of the scaling exponents describing the treeing process and their relationships to exponents 
related to the self-affine nature of the aggregate are reviewed. In section III a brief descrip- 
tion of both the RDSR and BD models and the numerical simulation technique is provided. 
The results are presented and discussed in section IV while the conclusions are stated in 
section V. 

II. Treeing and self-ctfRnity. Scaling relationships. 

Let us consider an aggregate that is growing above a d-dimensional substrate. Such 
aggregate is formed by trees that compete with each other leading to the entire pattern. The 
structural properties of both the trees and the resulting entire aggregate are determined by 
the growth mechanism. In these context, Racz and Vicsek^° have shown that for self-similar 
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fractals the tree size distribution that results from the competitive growth process is related 
to the structure of both the entire aggregate and the individual trees. Subsequently, this 
concept was extended to the case of self-affine objects^^. 

In the following paragraphs we will briefly outline the well known^°'^^ relevant definitions 
and the scaling relationship in order to establish the framework for the subsequent numerical 
study. Pointing our attention to dead trees of size s {s is the number of particles belonging 
to the tree) one has that both the rms height (hg) and the rms width (wg) of the trees obey 
simple power laws given by 

hscxs''^ (4) 

and 

Wg oc s^-^, (5) 

where and i'± are the correlation length exponents parallel and perpendicular to the grow- 
ing direction of the aggregate^^. These exponents may be different for self-affine (anisotropic) 
aggregates, while for self-similar objects i'\\=i'±. 

Assuming that Nt is the total number of particles of the aggregate, the average particle 
number N per unit area of the substrate is given by 

N = Nt/L'^, (6) 

where particles grow on a d-dimensional hypercubic substrate of size L in a D-dimensional 
space. 

Defining the number of trees with s particles Ns{N) in the whole aggregate one has that 

7Vt = J]sA^.(A^), (7) 

s 

and the cluster size distribution ngCN) that gives the probability of having trees of size s per 
unit area of the substrate is given by 

ng{N)^Ng{N)/L^. (8) 
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Replacing Eqs. (6) and (8) in (7) gives 



,10,11 



(9) 



s 



During the competition between trees along the evolution of the aggregate it may occur 
that the existence of large neighboring trees may inhibit the growing of smaller ones. This 
competing process ultimately leads to the death of some trees that become "frozen" within 
the underlying aggregate. These prevailing large trees continue the competition within 
more distant trees in a dynamic process. Since this situation takes place on all scales, it is 
reasonable to expect that the cluster size distribution should exhibit a power-law behavior 
so that 



where r is an exponent and f{y) is a scaling (cut off) function so that their asymptotic 
behavior is given by /(y) f=:i 1 for y <^ 1 and /(y) f^i for y » 1. 
Substituting (10) in (9) and after some algebra one can obtain^'^ 



Furthermore, recalling that the competition among trees actually is a dynamic process 
evolving in time {t oc Nt), one has that the survival time distribution of the trees (n^) may 
also obey a simple power-law behavior given by 



where 7 is an exponent. Therefore, one can establish a relationship between the size and 
the survival time t of the dead trees so that 



ns{N)r.s-^f{syN), 



(10) 



(7 = 2 - T. 



(11) 



rit oc t 



(12) 



s ^ t 



l/x 



(13) 



where the exponent x is given by 




(14) 
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Also, it is reasonable to expect that for the trees hs oc t, so inserting this relation in 
equation (13) and comparing it with equation (4) gives 

x = z/||. (15) 

Let us now establish a relationship between the growing trees and the properties of the 

aggregate. Assuming that the aggregate is self-similar (even compact) with fractal dimension 
D, its rms height h is expected to scale as 

hr^N''\ (16) 

for /i << L , where ui is an exponent^^. 

However, for length scales » h the aggregate is uniform in the remaining d lateral 
directions. This fact, in connection to equations (6) and (16) gives^^ 

D = i/-^ + d. (17) 

The rms height h of the aggregate can be rewritten as 

''' = ]^E"? = ]^E(E''?ww. (18) 

i s ies 

where hi is the distance of the i-th particle from the substrate. Taking into account that 
the average rms height of a tree of size s is given by 

hl = -^j:hl (19) 

equation (18) can be rewritten as 

h'-l^J2'h"ns{N), (20) 

s 

where equation (8) has been used. 

Replacing equations (4), (10) and (16) in (20) and performing some algebra one obtains 
the scahng relationship^^ 

2ui = -1 + (2 + 2z/|| - r)/a. (21) 



Now, using equation (11) and (17) it follows that 

a^v\\/n^v\\{D-d), (22) 

and substituting equation (22) in equation (11) one has^^ 

T ^2-v\\{D -d). (23) 

The scahng relationship (23) links the exponent r related to the tree size distribution 
to the correlation length exponent of the aggregate z/y . It is interesting to notice that is 
absent in equation (23) due to the fact that during the competitive growth process larger 
trees prevent smaller ones from growing. Scale invariance in the size distribution of the trees 
(equation (10)) is a result of this competitive process that is mainly governed by the tree 
hight (equation (4)). 

A typical example for the application of equation (23) is the growth of compact aggregates 
where D — d-\-l. So, equation (23) becomes 

T = 2-i.||. (24) 

Also, using equations (14), (15) and (23) one obtains 

7 = 1M- (25) 

There is no reason to expect that a compact aggregate may result from the addition of 
compact trees. In fact, the volume of trees of size s {vg) scales as 

- /i,^i;f-^ (26) 

and we define the following relationship between the volume and the number of particles 
~ s'^, where tt is an exponent, so one has 

v\\ + {D- = TT, (27) 

where for compact trees one has tt — 1 while for non-compact trees one has tt > 1. 
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Identifying the parallel correlation length with w{s) and the time with h{s), one has^^ 

z = (28) 

Using equations (14), (15), (23), (27) and (28) for trees in a compact aggregate {D — 
d + 1), we can obtain the following relationships 

u\\^Tiz/{z + D -1), (29) 

u±^7r/(z + D-l), (30) 

t = {{2-7t)z + 2{D-1))/{z + D-1), (31) 

and 

-f={z + D- 1)/{ttz) (32) 

All these relationships establish links among exponents that characterize the self-affine 
nature of the interface of the aggregate {z — a/13) and those corresponding to the description 
of the treeing process occurring in the bulk (i/||,i/x,T,7,7r). 

III. The RDSR and BD models. Simulation methods and mesoscopic equa- 
tions. 

In the random deposition with surface relaxation (RDSR) model a particle is released 
from a random position above the surface and falls vertically until it reaches the top of the 
selected column. Of course, such particle is initially located at a distance larger than the 
maximum height of the interface. The deposited particle is allowed to relax to a nearest 
neighbor column if the height of the neighboring column is lower than the one corresponding 
to the selected column. Further details on RDSR aggregates can be found in references^'*. 
Concerning the dynamics of tree formation, it is worth remembering the procedure used to 
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define a tree. In the case of the RDSR model, a newly deposited particle is assumed to 
belong to the tree corresponding to the impingement site, before its eventual relaxation. 

The ballistic deposition (BD) model is rather simple to describe: a particle is released 
from a random position above the interface of length L. Subsequently, the particle follows 
a straight vertical trajectory until it reaches the interface, whereupon it sticks. In contrast 
to the RDSR model, in the BD model no further relaxation of the particle is considered. 
Snapshot configurations of BD aggregates, and further details on the deposition rules can 
be found in references''^. Let us also remember that, in this case, trees are formed assuming 
that any newly deposited particle belongs to the same tree as that of the nearest neighbor 
particle where it is attached. If the deposited particle has more than one nearest neighbor 
belonging to different trees, one of them is selected at random and the particle is incorporated 
into that tree. 

RDSR and BD aggregates are grown in the direction perpendicular to d-dimensional 
substrate, i.e. in (d + l)-dimensions, using samples of different sizes (L*^) with 1 < d < 5. 
Simulation results are averaged over many different runs, depending on L and d. A Monte 
Carlo time step (mcs) involves the deposition of L'' particles. 

The interface of these aggregates is defined as the set of particles that are placed at the 
highest position of each column. So, the mean height of the interface (equation(l)) and the 
interface width (equation(2)) can be calculated and both are measured in lattice units (LU). 
It is well known that these aggregates are self-affine and the width of the growing interface 
w{L,t) scales as equation (3)"^. 

In contrast to the microscopic details of the growing mechanisms of the interface in both 
models, continuous equations focus on the macroscopic aspects of the roughness. Essentially, 
the aim is to follow the evolution of the coarse-grained height function h{x, t) using a well- 
established phenomenological approach that takes into account all the relevant processes 
that survive at a coarse-grained level. 

This procedure normally leads to stochastic nonlinear partial differential equations that 
may be written as foUows'^'^'^'^^'^^'^^ 
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= G,{/i(x, t)} + F + r?(x, t), (33) 

where the index i symbohcally denotes different processes, Gj{/i(x, t)} is a local functional 
that contains the various surface relaxation phenomena and only depends on the spatial 
derivatives of /i(x, t) since the growth process is determined by the local properties of the 
surface. Also, F denotes the mean deposition rate and r/(x, t) is the deposition noise that 
determines the fluctuations of the incoming flux around its mean value F. It is usually 
assumed that the noise is spatially and temporally uncorrelated, so fluctuations are given 
by a Gaussian white noise 

(r/(x,t))=0, (34) 

and 

(77(x, t)r]{^\ f)) = 2C5'^{x - x')5{t - t'), (35) 

where the brackets denote ensemble averaging, C is the strength of the fluctuations, and d 
the spatial dimension of the surface. 

The RDSR model can be described by the Edwards- Wilkinson equation given by^'^'^^'^^ 

= F + i.„V^/i(x, t) + 77(x, t), (36) 

where Uo plays the role of an effective surface tension, since the i/oV^/i(x, t) term tends to 
smooth the interface. Equation (36) can be solved in Fourier space and the following values 
of the exponents are obtained: z = 2, a = {2 — d)/2 and /3 = (2 — c?)/4. So, there exists an 
upper critical dimension dc — 2 above which one has a — Q and /3 = 0. 

The BD model can be described by the Kardar-Parisi-Zhang equation given by^ 

= F + vV'h + ^{Vhr + ry(x, t), (37) 

where A plays the role of an effective lateral growth. In this equation, the exponents satisfy 
the scaling relation 
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z + q; = 2, (38) 

as a consequence of Galilean invariance in the related Burgers equation^^. In the usual field 
theory approach, a perturbation theory is defined with respect to the nonlinear term (A). 
The corresponding renormalization group reveals^^ that the physics of the KPZ is related 
to a strong coupling fixed point, which is inaccessible by perturbative methods. Except for 
d = 1, where, thanks to a fiuctuation-dissipation theorem, an analytic solution is possible. 
In this case one has a = 1/2, (3 = 1/3, so that z = 3/2. Extensive numerical simulations 
have been carried out for restricted sohd on sohd models^^. These are discrete models that 
belong to the KPZ universahty class^° and these studies show a gradual decrease in the value 
of a, without any evidence for an upper critical dimension. In the same direction a recent 
study performed using a nonperturbative renormalization of the KPZ equation suggests a 
gradual decrease in the value of a^^. On the other hand, several analytical approaches'^ 
suggest an upper critical dimension (ic = 4 above which one should have a = 0. 

IV. Results and discussion. 

Figure 2 shows the tree distribution functions corresponding to the RDRS model in 
dimension 1 < d < 5. The tree size distribution (figure 2 (a)) for d = 1 exhibits a power law 
according to Eq.(lO) with exponent r = 1.35 ± 0.03. For higher dimensions Ug also exhibits 
power law behavior but the slope r = 1.50±0.01 becomes independent of the dimensionality. 
A similar behavior is shown by the survival time distribution (figure 2(b)) that has slopes 
7 = 1.54 ± 0.03 for (i = 1 and 7 = 2.00 ± 0.03 for d > 2, respectively. The measured 
exponents are listed in Table I. 

The distribution functions of the BD model also exhibit power law behavior (figure (3)). 
However, in contrast to the case of RDSR, for both Ug (figure 3(a)) and rit (figure 3(b)) the 
slopes depend on the dimensionality. The measured exponents are listed in the 5*^* and 9*'* 
columns of Table II. 
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The behavior of the rms height of the trees as a function of the tree size s is shown in 
figures 4(a) and 4(b) for RDSR and BD, respectively. Again, nice power laws are observed 
in all cases. For the RDSR model a clear change in the slope is observed between d = 1 and 
d >2, where for the higher dimensions the slopes are independent of the dimensionality (see 
S^'' column of Table I). For the BD model the slopes are different in all dimensions (see 3'"'' 
column of Table II) . 

Results shown in figure 5 for the rms width of the trees {wg) versus s are also fully 
consistent with the previous results and are hsted in the 2"'^ column of Tables I and II. 

Based on our numerical estimation of the exponents r, 7, i'± and i^w we are in a condition 
to check the scaling relationships discussed in section II that are summarized by equations 
(24), (27), (28), (31), (32) and (38), where the latter only holds for the BD model. 

The results obtained for the RDSR model (Table I) show that the exponents i'±, r 
and 7 are independent of the dimensionality of the surface for d > 2, as expected from figures 
1-4. These results are consistent with the fact that d — 2 is the upper critical dimension of 
the EW universality class. 

Inserting the exact value 2; = 2 for the EW universality class in equation (31) the 
"theoretical" estimation of r can be obtained (7*'* column of Table I), which is in excellent 
agreement with the numerical data. Also, inserting the measured values of i/y in equation 
(24) (see the 6*^* column of Table I) excellent agreement with the measured exponents is 
obtained. The quality of the achieved agreement is also excellent when comparing theoretical 
and estimated (equation (25)) values of 7, which is shown in the 10*^ and 9*^' columns of 
Table I, respectively. 

For the case of the RDSR model it is interesting to notice that the exponent tt (4*'* column 
of Table 1) departs from the trend exhibited by the other exponents. In fact, while tt ~ 1 
for d < 2 , such an exponent increases monotonically with the dimensionality for 0? > 2. 
Roughly, an average increment of 1/4 in tt for each additional dimension is estimated. The 
value of TT = 1 corresponds to the development of compact trees (sec equations (26), (27)), 
so it is concluded that for 0? > 2 the trees leading to the formation of RDSR aggregates 
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are non-compact objects. Therefore, the branches of the trees become interweaved forming 
complex patterns. 

Our results for the BD model are listed in Table II. In this case the measured exponents 
{i'±, r and 7 ) depend on the dimensionality, as expected from figures 2 — 5. The 
measured values of r (5*'* column of Table II) are in agreement with the estimation obtained 
using equation (24) (6*^^ column of Table II) for d < 4. For d = 5 this estimation presents a 
small deviation due to the large error involved in the evaluation of the exponent. However, a 
slight increment of r with the dimensionality is found in both cases (the measured exponent 
and equation (24)). This result is in contrast to the conjecture stating that d = 4 may be 
the upper critical dimension of the KPZ universality class^^. Furthermore, using equations 
(31), (38) and the values of a reported in reference [25] two additional comparisons can be 
made: inserting the values of z obtained using renormalization group calculations (numerical 
simulation) in equation (31) gives the values of r listed in the 7*'* (8*^) column of Table II. 
In both cases, excellent agreement with the measured exponents is found. 

We have to recognize large errors in our evaluation of 7 for the BD model (see figure 
2(b)) that prevent reliable estimations for d> A. However, comparisons with results obtained 
using equation (25) (lO^'^column of Table II), renormalization group calculation (11*'* column 
of Table II) and numerical evaluation (12*'' column of Table II), exhibit the trend already 
observed by our direct measurement of 7. 

It should also be noticed that tt ~ 1 (within error bars) pointing out that BD aggregates 
are built up by compact trees at least for 0? < 4, in contrast with our previous results ob- 
tained for RDSR aggregates for d > 2. 

CONCLUSIONS. 

Extensive numerical evidence is provided showing that the bulk of two typical growth 
models, belonging to both the KPZ and EW universality classes, can be rationalized in 
terms of a treeing process. Scaling relationships, linking exponents relevant to the treeing 
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process with standard dynamic exponents are reviewed and tested numerically in dimensions 
1 < d < 5. So, our main conclusions are summarized in Tables I and II where all the 
measured exponents are listed and compared to theoretical predictions. It is concluded that 
BD trees are compact objects (at least up to d = 4) while RDSR trees are non-compact 
objects for > 2 leading to interweaved complex structures. Based on the validity of 
equations (29-32) and since above the upper critical dimension the exponents i'±, r, and 
7 remain constant for self-similar aggregates, it can be stated that above the upper critical 
dimension the change in the dimensionality of the aggregate implies that tt must also change. 
Therefore, we conjectured that all kinds of self-affine aggregates above their upper critical 
dimension can be rationalized in term of the superposition of non-compact trees. 

Summing up, it is shown that the evolution of the treeing statistics is a suitable method 
for the characterization of the growing aggregates with promising applications in experimen- 
tal situations. In fact, it is well known that the electrochemical deposition method in thin 
cells provides a rather simple experimental setup for the observation of a wide variety of 
self-affine growing patterns^^'^^. In this case the dynamics of growing trees leading to frozen 
structures can easily be observed and the evaluation of the tree size distribution can be ob- 
tained straightforwardly after proper digitalization of the images of the aggregates. Another 
field of application is the study of poly crystalline thin film growth by scanning tunneling 
microscopy^^'^^ and atomic force microscopy^^. In these cases the direct evaluation of the 
number of crystallites as a function of the deposit height accounts for the number of frozen 
crystallites (i.e. the trees of the growing aggregate) and the crystallite size distribution can 
be evaluated. 
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Table I: List of exponents u±, i/y, r (5*'* column) and 7 (8*'* column) measured for the 
RDSR model for dimensions 1 < d < 5. The values of theoretical estimations of r and 7 
are obtained taking z = 2 in Equation (31) (7*'* column) and Equation (32) (lO*'^ column), 
respectively. For these estimations tt = 1 for d < 2 and tt = 5/4, 3/2, 7/4 for d — 3,4,5, 
have been used, respectively. The values of tt obtained using Equation (27) are shown in 
the 4*'* column. The values of r and 7 obtained using Equations (24) and (25) are shown in 
the 6*'* and 9*'^ columns, respectively. 



D 




Z/|| 


vr 


r 


r(Eq.24) 


r(t) 


7 


7(Eq.25) 


7(t) 


1+1 


0.35(2) 


0.63(1) 


0.98(2) 


1.35(3) 


1.37(1) 


4/3 


1.54(3) 


1.59(3) 


3/2 


2+1 


0.27(2) 


0.50(1) 


1.02(4) 


1.49(1) 


1.50(1) 


3/2 


1.98(5) 


2.00(4) 


2 


3+1 


0.26(2) 


0.50(1) 


1.28(6) 


1.50(3) 


1.50(1) 


3/2 


2.00(2) 


2.00(4) 


2 


4+1 


0.25(2) 


0.50(1) 


1.50(8) 


1.50(1) 


1.50(1) 


3/2 


1.99(2) 


2.00(4) 


2 


5+1 


0.25(1) 


0.49(1) 


1.74(4) 


1.50(2) 


1.51(1) 


3/2 


1.99(4) 


2.04(4) 


2 
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Table II: List of exponents v^, r (5*'* column) and 7 (9*'* column) measured for the 
BD model for dimensions 1 < c? < 5. The estimations performed using Equations (24) and 
(25) are shown in the and 10*^ columns, respectively. The exponents obtained using 
Equations (31), (32) and (38), tt = 1 and a values obtained by renormalization group 
calculation (RG) and numerical simulations (NS) reported in reference^^ are shown in the 
7*^, S*'^, 11*'', 12*'* columns. The exponent tt is evaluated using equation (27) (4*^* column). 



D 




Z/|| 


TT 


r 


r(Eq.24) 


r(RG) 


r(NS) 


7 


7(Eq.25) 


7(RG) 


7(NS) 


1+1 


0.40(1) 


0.60(1) 


1.00(2) 


1.40(1) 


1.40(1) 


1.40 


1.40 


1.70(2) 


1.66(3) 


1.66 


1.66 


2+1 


0.29(1) 


0.45(1) 


1.03(2) 


1.57(1) 


1.55(1) 


1.54 


1.55 


2.3(1) 


2.22(5) 


2.22 


2.24 


3+1 


0.23(2) 


0.38(1) 


1.07(6) 


1.65(1) 


1.62(1) 


1.63 


1.64 


3.2(2) 


2.63(7) 


2.75 


2.77 


4+1 


0.21(4) 


0.34(2) 


1.18(16) 


1.70(2) 


1.66(2) 


1.69 
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FIG. 1. Schematic view of the deposition of particles forming trees. In (a) the particles are 
deposited with the rules of ballistic deposition. Particle A belongs to tree 1 or 2 with the same 
probability, particle B belongs to tree 2 and particles C and D belong to tree 4. The growth of 
tree 3 is stopped after particle C deposition and it becomes a frozen tree. In (b) the particles are 
deposited with the rules of random deposition with surface relaxation. Particle A, which belongs 
to tree 2, can either relax to the right or to the left with the same probability. Notice that in the 
former, particle A partially shadows tree 1, while in the latter it simply becomes attached to tree 
2. Particle B belongs to tree 2 but it relaxes on top of tree 3 causing partial shadowing. Particle 
C belongs to tree 4 and after relaxation on top of tree 3 causes the growth of that tree to stop, so 
that it becomes a frozen tree. Finally, particle D belongs to tree 4. 
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FIG. 2. Log-log plots of the tree distribution functions corresponding to the RDRS model in 
dimensions 1 < d < 5. (a) The tree size distribution, where the size s is given by the number of 
particles of each tree, (b) The survival time distribution, where the time is measured in mcs. 
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FIG. 3. Log-log plots of the tree distribution functions corresponding to the BD model in 
dimensions 1 < d < 5. (a) The tree size distribution where the size s is given by the number of 
particles of each tree, (b) The survival time distribution, where the time is measured in mcs. 
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FIG. 4. (a) and (b) Log-log plots of rms height of the trees as a function of the tree size s for 
the RDSR model and the BD model, respectively. The rms height is measured in LU and the size 
s is given by the number of particles of each tree. 
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FIG. 5. (a) Log-log plots of rms width of the trees as a function of the tree size s for the RDSR 
model, (b) Log-log plots of rms width of the trees as a function of the tree size s for the BD model. 
The rms width is measured in LU and the size s is given by the number of particles of each tree. 
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